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Abstract 

In their original study of conformal gravity, a candidate alternate gravitational theory, Mannheim 

o 

CN ■ and Kazanas showed that in any empty vacuum region exterior to a localized static spherically 

H , 

symmetric gravitational source, the geometry would reduce to the standard attractive gravity 



Schwarzschild geometry on solar system distance scales. In a recent paper Flanagan has argued 
that this would not be the case if the source has associated with it a macroscopic scalar held which 
makes a non-zero contribution to the energy-momentum tensor in the otherwise empty exterior 
region. In this paper we examine Flanagan's analysis and show that even with such long range 



> 

m 
o 

scalar fields, the standard Schwarzschild phenomenology is still recovered. 

o 
7' 

u 

bX). 

> 



Electronic address: philip.mannheim@uconn.edu 



I. INTRODUCTION 



As a theory of gravity, conformal gravity is very appealing since it not only possesses 
the full metric structure associated with the standard Einstein gravitational theory, it in 
addition also possesses the local conformal symmetry which is characteristic of theories of 
elementary particles in which particle masses are to be generated entirely via dynamically 
generated symmetry breaking phase transitions in the vacuum. With the action of the 
conformal theory being required to be invariant under local conformal changes of the metric 
of the form g^(x) — > e 2a ^ g fll/ (x), in the conformal theory the purely gravitational sector of 
the action is uniquely prescribed to be of the form 



Iw = -oc 9 j d 4 x(-g) l / 2 C x ^ K C x ^ K 

= -2a g J d 4 x(- g y/ 2 \r„ k r» k - 1 -{r 



a I 



(1) 



where C\ I1VK is the conformal Weyl tensor and a g is a necessarily dimensionless gravitational 
coupling constant. Similarly, with the matter action equally needing to be conformally 
coupled at the level of the Lagrangian, the prototypical conformal matter action consisting 
of a fermion field ip(x) with spin connection T^x) and a scalar field S(x) is uniquely specified 
to be of the form 
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J d^-g) 1 ' 2 [-S^S. 4l - ^S 2 R^ + \S A + i^Y(x)[d, + F,(x)]i> - hSW 



(2) 



where the coefficient of the S 2 R' 1 ll term is uniquely required to be equal to the indicated 
negative factor of —1/12 and where h and A are dimensionless coupling constants. On 
the theoretical side, giving the scalar field S(x) a non-zero vacuum expectation value will 
spontaneously break the conformal symmetry and give the fermion a mass (as needed for 
particle physics), while on the phenomenological side the conformal theory has been found 
capable of readily solving both the dark matter and dark energy problems which currently 
challenge the standard theory (see e.g. [l] where full bibliography and background are given). 

For the above Iw + Im action functional variation with respect to the matter fields yields 
the equations of motion 

zr(x)[d, + T,(x)]^-hS^ = 0, (3) 



and 



S ^ + l S R%-AXS 3 + h^ = } (4) 



while functional variation with respect to the metric yields 

AagW^ = T^, (5) 

where W^ v is given by 

-lg^(R a a y% + \{R\)^ V + \r\R" v - V(#\J 2 , (6) 



and T^ v is given by 

= i^y*(x)[9" + T v (x)U + -S^S> V - -g^S' a S. a - -SS> KU + -g^SS' a a 

3 6 ' 3 3 

As 2 (V - \g» v R a ^ - g» u \S\ (7) 

For the purposes of studying localized static spherically symmetric matter sources we can 
represent the fermionic part T^ u = i%p r y^ \x)[d v + T V '(x)}ip of the full energy-momentum tensor 
T^ u as the perfect fluid 

f»» = ~[(p + p)U> l U v +pg^}. (8) 
With the trace of this T^ v being given by 

n = \ [3p - p] = i^(x)[d, + = hSW, (9) 

we find that Eq. (@J can be rewritten as 

l -S 2 R^ = - c [p - 3p] + 4AS 4 - SS-%, (10) 

a form that will prove central in the following. 

n 

Given the form of W^ v in Eq. (EJ), it can immediately be shown [2j that in any p(r > 
R) = 0, p(r > R) = 0, S(r > R) = empty, static, spherically symmetric exterior r > R 
region where the full T^ v vanishes, the R^ v = Schwarzschild solution is indeed an exterior 
empty vacuum solution (W^ = 0) to the theory; with its matching to the interior region 
at the surface r = R showing [3] that the coefficient 2(3 of the — 1/r term in the familiar 
Schwarzschild solution — goo — VflVr = 1 — 2(3 /r is indeed the positive one needed for 
gravitational attraction. If however, we instead allow the scalar field to be non-zero in the 
r > R region where the localized p(r) and p(r) are to vanish, on making the specific local 
conformal transformation S(r) — > e~ a(j '^S(r) which brings the scalar field to a constant value 



So in such r > R exterior regions, and on then, as per the analysis followed in 4|, dropping 
the ensuing —q^XSq term from Eq. ([7]), the full T^ u is then reduced to 

T^(r) = \{{P + vWU v + pgn - ^S 2 (rt - l -g^R a ^j , 

T^(r > R) = ~S 2 (rT - \g^R a ^j , (11) 

to now not necessarily vanish outside the source. However, before proceeding to analyze the 
implications of Eq. (ITT]) , we note that the dropping of the —g^XS^ term while retaining 
the — (jSq/6) (^R^ — ^g^R ^ term is not actually consistent with the notion of symmetry 
breaking in curved spacetime, since the very existence of a constant scalar field which is 
non-zero throughout the spacetime is a global, long range order effect and not a local one, 
an effect which thus couples to the global geometry. Indeed, it is the very non-vanishing 
of the —g^XSQ term in which is needed to support the non-zero So in the first place, 
with the AS 4 - (S 2 /12)R a a term in the action of Eq. Q2) serving as a double-well potential 
in spacetimes (such as for instance de Sitter) in which the global curvature is non-zero. 
The —(Sq/6) {^R^ v — -g^R 01 ^ term in T^ u is thus just as much a part of the cosmological 
background as the — g^XS^ term. In fact, as we shall discuss in detail below, the effect 
of the introduction into the cosmological background of a local mass source is not actually 
associated with the constancy of the scalar field, but rather with the departure from con- 
stancy which is brought about by the very presence of the localized mass source. And it will 
actually be this spatially dependent departure from constant field and not the — (Sl/12)R a a 
term itself which will be responsible for the gravity produced by a localized source. 

Nonetheless, even with this proviso, it is still of interest to determine what particular local 
gravity would be associated with the truncated energy- momentum tensor given in Eq. ( TTTT) . 
and to this end we note that in the exterior region an R^ v = geometry will nonetheless 
actually cause the T^ v [r > R) of Eq. (TTTj) to vanish, with an R^ u = geometry still being 
an exact = 0, T^ v = solution to the theory in the r > R region. However, with Eq. 
(fit)]) reducing to 

^R% = - c [p-3p], (12) 

and with the standard Einstein equations for a standard pure T^ v perfect fluid source and 
no macroscopic scalar field, viz. 

(V - lsrx*a) = - c i(p+ pWW + P9" u ] , (13) 



8ttG 



entailing that 



-—BP 

8ttG " 



1 



[3p - p] , 



(14) 



the flip in sign between Eqs. (12) and (14) (as occasioned by the —1/12 factor in Eq. (j2J)) 
would suggest that the conformal gravity B^ = Schwarzschild solution would correspond 
to repulsive rather than attractive gravity. Precisely such a concern has been raised by 
Flanagan J4] and Van Acoleyen [5], to thus call into question the viability of the conformal 



theory. However, an explicit non-perturbative numerical analysis of Eqs. 



n 



J5J by Wood 



and Moreau [6J has found that even with a macroscopic scalar field, solar system gravity 



is still attractive. It is the purpose o 
particular to examine the analysis of 



)f this paper to reconcile these differing claims, and in 
4| in detail. However, in order to do this we need first 
to recall how it was that Mannheim and Kazanas were able to recover standard attractive 
gravity in an empty exterior region in the first place. 



II. EMPTY EXTERIOR REGION SOLUTION 



To deal with the case of a static spherically symmetric geometry, rather than use a metric 
in the standard form 

ds 2 = -b(p)c 2 dt 2 + a(p)dp 2 + p 2 dVt 2 , (15) 
in the conformal theory it is convenient to make the general coordinate transformation 



P = p(r), B(r) = ^ A(r) = r*a(r)p*(r) 



(16) 



p 2 (r) ' 1 p 2 (r) 
with an initially arbitrary function p(r), as this brings the metric of Eq. (fl5l) to the form 

p 2 (r) 



ds 2 



On now choosing p(r) according to 



B(r)c 2 dt 2 + A(r)dr 2 + r 2 dVl 2 



dr 



(17) 



(18) 



p(r) J r^G^r)^?")] 1 / 2 ' 
the function A(r) is made equal to 1/B(r), with the metric of Eq. (fT8|) then being found to 
take the form 

p 2 {r) 



ds' 



- B (r)c 2 dt 2 + ^- + r 2 dn 2 
K ' B r) 



(19) 



As such, the above sequence of purely kinematic coordinate transformations trades the two 
metric coefficients a(p) and b(p) for two other coefficients B{r) and p(r), but does so in 



way in which one of them, p(r), appears purely as an overall multiplier of the metric. Since 
the theory is conformal invariant, both of the functions W^ v and T^ u which appear in Eq. 
§5§ transform under the conformal transformation g^ v — > e 2a( - x ^g^ as — > e~ 6a( - x ^W fJ " u , 

can thus be scaled out of the theory, with the full 
content of the conformal theory being contained in the metric 

dv 2 

ds 2 = -B{r)c 2 dt 2 + — + r 2 dQ 2 , (20) 
a metric which now contains just a single unknown metric coefficient B(r). 

fin 

The great utility of the metric of Eq. (f20l) is that for it yj the combination W° — W r r 
evaluates exactly and without any approximation whatsoever to the extraordinarily compact 
form 

±(W° -W r r )=B""+—. (21) 

(The primes denote derivatives with respect to r.) Consequently, on recognizing that 
B"" + 4B'"/r is the radial piece of the fourth order Laplacian V 4 -B, and recalling that 
the tracelessness and covariant conservation of W^ u \W° + W T T + 2W% = 0, (B'/2B — 
l/r)(W° — W r r ) — (d/dr + 4/r)W r r = 0] entail that it only has one independent component 
in the static spherically symmetric case, we find that the equation of motion of Eq. (JS]) can 
be written very compactly as 

V 4 5(r) = /(r) (22) 
where the source function f(r) is defined via 



ia s B V u 'I 4a s 



3(P + P) vSS „_ 2S a 



cB ... (23) 

as evaluated here with the metric given in Eq. ( l20l) and the explicit form for T^ v given in 
Eq. ©. 

As a differential equation, Eq. ( 122]) admits of the general solution 

f rr 1 rr 1 roo y*2 poo 

B(r) = — / dr'r' 2 f(r') / dr'r' 4 f(r') - - / dr'r ,3 f(r') / dr'r'f(r'), (24) 

2 jo 6r jo 2 J r 6 J r 

together with an additional cosmologically relevant w — kr 2 term which solves the homoge- 
neous V 4 -B(r) = 0. Thus, in cases in which the entire source f(r) vanishes in the r > R 
region (viz. the S(r > R) = 0, p(r > R) = 0, p(r > R) = case considered in the 
metric exterior to the source is given as 

28 

B(r > R) = 1 - -j- + 7r (25) 
6 



where we have introduced the coefficients 

7 = -~ f R dr'r' 2 f(r'), 2(3 = \ [* dr'r*f(r'). (26) 
2 Jo b Jo 

With the (3 and 7 coefficients being associated with different moments of the source they are 
in principle different, with the successful phenomenological fitting of conformal gravity to 
galactic rotation curve data |l( without the use of any dark matter leading to the numerical 
values (3 = 1.48 x 10 5 cm, 7 = 5.42 x 10 -41 cm -1 per unit solar mass of material. With such 
values, the linear potential modification to the Newtonian 1/r potential is thus negligible on 
solar system distance scales, with the standard solar system Schwarzschild phenomenology 
then following. For these particular values of (3 and 7 we note also that the ratio (P/'j) 1 ^ 2 
is equal to 0.52 x 10 23 cm, a ratio which is thus galactic in scale rather than stellar. Hence 
the function f(r) which appears in Eq. f)26p must be singular, since if it were uniform it 
would otherwise lead to (/3/7) 1//2 ~ R for a star of radius R. Thus unlike the second order 
Poisson equation which is not sensitive to the structure of the source, precisely because it is 
a higher derivative equation, the fourth order Poisson equation does explore the singularity 
structure of the source. To get a sense of what the sensitivity to singularities might be, we 
note that the illustrative (though not mandatory) choice of source 

5(r) 



f(r) = -2p 
leads to 



..2 



3q 






|"<5(r)l 


2 








Q, 


J = p. 





(27) 



2(3 = q, j = p. (28) 

With the p-dependent term only contributing to the second moment and the g-dependent 
term only contributing to the fourth, the source of Eq. ( 1271) shows that in the presence 
of singularities the second and fourth moments are logically independent, a point that will 
prove to be of relevance in the discussion below. (As a source we note that its g-dependent 
part can be written as the e — > limit of the quantity 6ge(9r 4 — 10r 2 e 2 — 3e 4 )/7r(r 2 + e 2 ) 5 , 
a quantity which (for q > 0) is positive in r > 0(e) and negative in r < 0(e), but which 
traps a singularity at r = in the e — limit, causing its net contribution to the second 
moment integral of the source to vanish. Such singularities while foreign to the conventional 
standard gravity wisdom are not forbidden by any known gravitational observation - even 
the cherished positivity of the energy density is not actually mandated by observation, as 
the positivity of a moment integral does not entail the positivity of its integrand.) With the 



sign of (3 being fixed by the sign of q and with the sign of q itself being fixed by the sign 
of ot g in Eq. (1231) . there thus will be a choice of sign for the coupling constant a g for which 
the coefficient (3 will be positive, with it thus being possible to unambiguously fix the sign 
of the coefficient of the 1/r term in the conformal theory. (In the standard theory this same 
outcome is achieved by an a priori choice for the sign of G.) Additionally we note that the 
sign of a g bears no relation to the sign of the coefficient of the S 2 R^^ term in the action of 
Eq. ([2]), a point which will also prove to be of relevance below. 

Proceeding next to the case where there is to be a non-vanishing S(r) in the r > R region, 
we note that with the scalar field transforming as S(x) — > e~ a ^S(x) under a conformal 
transformation, one should not expect that the transformation which will bring the metric to 
the form of Eq. f[2T)j) will be just the one which will completely cancel the spatial dependence 
of S(r) and bring it to a constant So- Consequently, if one takes the metric to be of the 
form given in Eq. (1201) . one should not expect the associated T^ u to contain no derivatives 
of S(r). Working with the metric in the form of Eq. ( 1201) thus obligates working with the 
full spatially dependent S(r) in Eq. (JTj). For the particular choice of metric of Eq. (1201) Eq. 
(TlOl) is found to take the form 



-S 2 
6 



+ 4AS 4 - S 



BS" + — + B'S' 



(29) 



which when taken in conjunction with Eq. (122]) and an equation of state for the perfect fluid 
thus defines the problem in the non-vanishing S(r) case. Unfortunately no exact analytic 
solution to this set of coupled equations is currently known, and for the moment one must 
resort to a numerical treatment. A numerical solution has been provided in 6| which finds 
that the standard attractive Schwarzschild phenomenology continues to hold in the solar 
system, with the wisdom of the solution of Eq. ( 1251) thus being maintained. (Numerically, 
the effect of the gradient of the scalar field on both the geometry and particle trajectories 
in the exterior region was found to be within current experimental bounds, while leading 
to a small departure from the standard Schwarzschild phenomenology which the authors of 
jfl suggest could be responsible for the anomalous acceleration of the Pioneer spacecraft - 
the radial gradient of any macroscopic scalar field associated with the sun would naturally 
point toward it.) Armed with the above analysis we turn now to the work of Flanagan. 
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III. LINEARIZATION IN ISOTROPIC COORDINATES 



With the metric associated the exterior empty vacuum solution of Eq. (|25|) being close 
to flat in the intermediate region where the radial distance r is neither too small nor too 
large, in this intermediate region one can thus consider solving for the case in which there is 
a constant exterior macroscopic scalar field by linearizing the theory around flat spacetime. 
This then is the approach of Flanagan, and since one would not expect that it would matter 
what choice of perturbed coordinate system one might choose to work in, Flanagan opted to 
work with isotropic coordinates rather than standard ones. However, as we shall see below, 
on working the problem through in linearized standard coordinates we will actually reach 
a conclusion quite different from the one that Flanagan reached in his linearized isotropic 
coordinate study. We shall thus describe both of the linearization calculations and shall 
then reconcile the difference between them. 

For an isotropic coordinate system with metric 



ds 2 



H(p)c 2 dt 2 + J(p) dp 2 + p 2 d£l 2 



(30) 



linearization is associated with the metric 



d s 2 = -[1 + h(p)]c 2 dt 2 + [l+j(p)\ dp 2 + p 2 dVt 2 



(31) 



where h(p) and j(p) are small. With the explicit form of the exact W^ u associated with the 
metric of Eq. (J30l) being given in [7], its linearization is found to lead to 

1, 



\ir-h""] + ^if'-h>"] = ^[j-h], 
3 Sp 3 



(32) 



3p 3p 

Similarly a linearization of the Einstein tensor G^ u = — (l/2)g flu R a Q in isotropic coor 
dinates yields 



j" + -J 
P 



-V 2 j 



-~\j' + ti], 



(33) 
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with the Ricci scalar being given by 



R\ = \2j" + h" 



P 



[2f + ti] = V 2 [2j + h] . 



(34) 



Thus in the non-relativistic perfect fluid limit where p(p) <C p(p), Eqs. ([5]), (TTTT) and (|T2|) 
reduce to 

An, n(n\ R 2 _ 

(35) 



4a.W = %V* b - - ft] = £<« - f V*J, 
3 c 6 



4rv 

4a„W» = — 9 - 



p p 2 p A 

^-R^ = ^[2j + h] pip) 



C2 



(36) 
(37) 



6^6 c 
In his paper Flanagan explored various classes of solution to this set of equations, but if we 
are to recover the linearized isotropic coordinate system Schwarzschild solution, viz. 



ds 2 = - 



' 2MG~ 




2MG 




c 2 dt 2 + 


1 + 


P 




P 



dp 2 + p 2 dVl^ 



(38) 



we should set j(p) = —h(p), with Eq. (37) then leading to the wrong sign for the coefficient 
of the 1/p term. This then is the concern raised by Flanagan. 

IV. LINEARIZATION IN STANDARD COORDINATES 



To assess the significance of this result it is instructive to repeat the analysis in standard 
coordinates. For the standard coordinate system with metric 



ds 2 = -B(r)c 2 dt 2 + A(r)dr 2 + r 2 dfl 2 , 
linearization is associated with the metric 

ds 2 = -[1 + b(r)]c 2 dt 2 + [1 + a{r)]dr 2 + r 2 dVL 2 . 
This time the linearization of the exact W^ v given in [7| is found to lead to 

W tl = 



1 


b"" 


4b'" 
+ + 


a'" 


3 




r 


r 


1 


~b'" 


+ 2b" 


2b' 


3 


r 





2a' 2a 



2a 



(39) 



(40) 



(41) 
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Similarly, a linearization of the Einstein tensor in standard coordinates yields 



a a 



G\. 

with the Ricci scalar being given by 



b' a 

b" V a' 
"2 ~ 27 + 2? 



(42) 



R?,, = b" + 



2b' 2a! 2a 

ry> iyi ly 2 



(43) 



Thus in the non-relativistic perfect fluid limit where p(r) <C p(r), Eqs. (131) , (iTTj) and (TI2|) 
reduce to 



Aa g W tt 



b"" 



Ab'" a'" a" 2a' 2a 
+ + — + -- — + — 

ry ry iy ^ ry<J iy rt 



c 6 



Aa g W rr 



4a, 



b'" 2b" 2b' a" 2a 

w ry*2i jy3 ty2 y4 



C2 C2 

6 6 



b „ 2b' 2a! 2a 

iy /y ty*2i 



C2 

' 6 
P 



b' a' 



(44) 

(45) 
(46) 



to thus define the problem in the standard coordinate system. 

Before attempting to look for solutions to Eqs. (|44l) - (|46l) it is useful to recall how the 
linearized standard coordinate Schwarzschild solution is obtained in the standard Einstein 
theory. Specifically one linearizes Eqs. ([IB"]) and (fill) using the metric of Eq. (1401 and 
ignores the perfect fluid pressure, to obtain 

^3 



87TG 



8ttG 



b" + 



2b' 2a! 



a a 
- + - 

ry :y ij 

2a 



P 
■ 

c 
P 

j 

c 



VIZ. 



8llG 

with general all r solution 

a(r) 
b(r) 



a a 
- + - 

iy y £• 



8ttG 



6" + 



2b' 



(47) 



(48) 



8ttG 



<rr Jo 
8ttG 



err Jo 



dr'r' 2 p(r'), 
dr'r' 2 p{r') 



8ttG 



dr'r'p(r'), 



(49) 
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and familiar r > R exterior solution 

2MG , , 2MG , , 

b(r>R) = — , a (r>i2) = — =— 50 

err err 

for a source with rest energy Me 2 = Air J R dr'r' 2 p{r'). 

With Eq. (j4*6l) only differing from the second equation in Eq. (j47H by the replacement of 
.Sq/6 by — c 3 /8nG, any attempt to solve Eqs. (jHJ) - (Pf6l) via an analog of Eq. (1491) would 
immediately lead to us a Schwarzschild geometry with repulsive gravity. Since we wish to 
recover the exterior attractive gravity Eq. (1501) in the conformal case, we must thus not 
recover the analog of the interior Eq. ( H9j) as well, i.e. in the conformal case we need to find 
some different r < R solution which will still match on to the exterior Eq. fl50|) . Noting now 
that Eq. (H7|) involves V 2 b but not V 2 a, we see that the a(r) sector of the standard gravity 
Eq. (f4"Tj) does not involve any V 2 (l/r) term with its associated — 4ir5(r)/r 2 singularity, with 
the 1/r term in a(r) not needing any delta function source to support it. Consequently, 
returning now to the conformal case where there is equally no V 2 a term in Eq. (T4"6T) . in the 
conformal case we can thus look for a solution in which a(r) has a 1/r form for all r, both 
exterior and interior. For all r we thus set 

a ( r ) = J; - er ( 51 ) 

where d and e are as yet undetermined constants, and where the er term that we also 
introduce will prove useful below. (With Eq. (|44p containing no V 4 a term, the er term is able 
to be present in a(r) at all r also, as it too will not generate any delta function singularity.) 
Quite remarkably, we find that for our candidate form for a(r), the contribution of its 2d/c 2 r 
term not only cancels identically in Eq. (T4*6l) . it also cancels identically in Eq. (j4"4"l) as well. 
In addition, the er term also drops out of the left-hand side of Eq. (jHJ) identically, though 
it does generate terms of the form eS^/r elsewhere in these same equations. However, we 
shall show below that on solar system distance scales such eS 2 /r terms are totally negligible 
compared to the p/c terms in Eqs. f|44|) and f|4*6l) . with their neglect then bringing Eqs. ( l4*4"j) 
and fj4*6|) to the form 

VW* V%=%. (52) 
4a g c cSq 

With the generic forms of the solutions to these Poisson equations being given as in Eqs. 
and (|49p . the exterior solution thus has to simultaneously be given as both 

Mr >R) = ^— / dr'r' 2 p{r') + —— / dr'r' A p{r ! ) (53) 
8a g c Jo 8a g cr Jo 
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b(r > R) = --J- [ R dr'r' 2 p{r'). (54) 



and 

hfr ^ fi^ — _ 

^ocr Jo 

With the rest energy of the source being given by the positive Mc 2 = 4tt J r dr'r' 2 p(r'), the 
solution of Eq. (154"|) immediately assures us that its associated Newtonian potential is indeed 
attractive. Compatibility of the Newtonian terms contained in the solutions of Eqs. (|53|) 
and (l54"j) requires that the involved parameters be related via 

6Mc M 4 c 

AnSl ~ ~32na g ' ^ ' 

(we use M4C 2 = Air dr'r' 4 p(r') to denote the fourth moment of the energy density), with 
the right-hand side of Eq. (|55|) indeed being able to be positive for an appropriate choice of 
the sign of a g . Compatibility of Eqs. f[53"j) and (15"41) additionally requires that the coefficient 
a g be large enough to make the linear potential term 3rMcj32 r na g be negligible on solar 
system distance scales where we are looking for a linearized solution in the first place (i.e. as 
we had noted previously, since the fourth order theory leads to both Newtonian and linearly 
rising potentials, we can only linearize in the region where r is neither too small nor too 
big). Finally, on inserting the exterior solution of Eq. ( 1531) for b(r) into Eq. ( 1451) . the one 
equation still remaining, we see that in this solution the left-hand side of Eq. (1451) will vanish 
identically in the exterior region if we take the strengths of the linear potential terms in the 
exterior a(r) and b(r) to be equal and opposite (viz. e = 3Mc/32ira g ), with the associated 
exterior vanishing of the right-hand side of Eq. (|45j) (viz. none other than the vanishing 
of G r r in the exterior region) then fixing the parameter 2d/c 2 in the expression for a(r) 
given in Eq. ( 15T1) to be of the form 2d/c 2 = 6Mc/47tSq, just as required for the attractive 
exterior Schwarzschild solution. (The contributions of the linear potential term in b(r) (viz. 
b[r) = er) to Eqs. fj45|) and fj46|) are also suppressed in regions where S^e/r is negligible.) 
It is thus the freedom to assign both the magnitude and the sign of a g which enables us 
to achieve compatibility of the solutions of Eqs. fl53|) and (1541 on the solar system distance 
scale region of interest to us. 

To show that there is at least one explicit choice for p(r) for which the needed compat- 
ibility between the solutions of Eqs. ( |53l) and (1541) can be made manifest, guided by the 
highly singular structure of the source given in Eq. (I2T1) we set 



, , Mc 2 5(r) 3vc 



2 



(56) 
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to obtain as the two respective forms for b(r > R) 

3Mcr 3v 

b(r > R) = (57) 

v ; 32na g Aa g r y ' 

and 

^ > fl » = -sap < 58 > 



In Eq. p8j) , and central to the analysis here, we note that the carefully chosen v -dependent 
term in the source given in Eq. (1561) makes no contribution at all to the second moment 
integral needed in Eq. (154"|) . with the second order Poisson equation simply being totally 
insensitive to its possible existence. The identification of a local effective gravitational 
constant Gloc according to 

" (59) 



Si 3c 3 

allows us to write Eq. (l58l) and the 1/r part of Eq. ( IBTl) as the attractive exterior metric 

Hr > R) = „(r > «) = ^ (60) 



Compatibility of this solution with that of Eq. (1571) requires that the w and a 9 parameters 
be constrained according to 

3Mc _ 2MG LOC _ 3t; 

2tt^ 2 c 2 4a fl ' 1 J 

a relationship which can be satisfied by making an appropriate choice for the magnitude and 
sign ofv/otg. (In a full non-perturbative treatment of Eqs. (jSj) and (JTUj) the radial dependence 
of the energy density and pressure of the fluid would be self-consistently determined as part 
of the solution. Here we seek only to show that in the linearized case there is a choice for 
p(r) for which one can obtain compatibility between the solutions of Eqs. (I5TI) and (1581) . ) 

As regards the numerical values of the various coefficients in our solution, for one solar 
mass of material we want the solution of Eq. (1571) to recover the B(r > R) = 1 — 2(3 /r + 
solution whose fitting to galactic rotation curves described earlier requires that the coeffi- 
cients take the numerical values 2(3 = 2.96 x 10 5 cm, 7 = 5.42 x 10 -41 cm -1 . For these 
particular values we find that if we take Gloc to be equal to the standard Newtonian G, the 
parameters in Eqs. (1571) and (|58|) are then given as Sq = 9.66 x 10 37 gm sec -1 , v = 1.29 x 10 88 
gm cm 3 sec -1 , a g = 3.29 x 10 82 gm cm 2 sec -1 (i.e. a g has the dimension of action as re- 
quired to make it dimensionless in natural units); and with these numbers we find that 
the 3Mcr /32iia g and eS^/r terms are indeed completely negligible on solar system distance 
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scales where we are seeking a linearization of the theory (viz. 3Mcr /32na g <C 3Mc/2irS 2 r 
and eS 2 /r <C p/c ~ Mc/R 3 ). Our analysis thus shows that even with a constant macro- 
scopic scalar field in the region exterior to the sun, it is still possible for conformal gravity 
to be compatible with solar system Schwarzschild phenomenology, with it being central to 
our analysis that the source possess much deeper singularities than the sources that are 
ordinarily considered in the standard theory. 



V. RECONCILIATION OF THE CALCULATIONS 

Since we have constructed an explicit solution for the linearized metric coefficients b(r) 
and a(r) using the standard coordinate system, it must be the case that on making a coordi- 
nate transform to isotropic coordinates the resulting metric coefficients h(p) and j(p) must 
provide a solution to the problem in the isotropic coordinate system. And yet inspection 
of Eq. ([37]) would suggest otherwise. To reconcile this seeming contradiction we note that 
while the b(r) and a(r) coefficients respectively transform into h(p) and j(p), because the 
transformation between the full non- linearized metrics of Eqs. ( 1301) and (13T?|) is of the form 

r = pJ^(p), ^ 1/2 W = 1 + f^> B(r) = H(p), (62) 

in this transformation A(r) is related not to J(p) but to a derivative of it instead. Thus, 
what had initially been a first derivative function of a(r) in Eq. (I46p becomes a second 
derivative function of j(p) in Eq. (1371) . Specifically, with the linearized transform itself 
being of the form 

r^p+^-, - a (r)^pj'(j>), b{r) - h(j>), (63) 

the -2a /r 2 - 2a! jr combination which appears in the expression for the Ricci scalar given 
in Eq. (]46l) thus transforms into the 2V 2 j(p) term which appears in the expression for the 
Ricci scalar given in Eq. (13T1) . with Eq. (1631 explicitly leading to 



2a! 2a 



f + 2f 



2V 2 j (64) 

in lowest order. Now the 1/r part of our chosen standard coordinate system solution for a(r) 
as given by a(r) = 2d/c 2 r has the property that it causes the quantity —2a/r 2 — 2a! jr to 
vanish identically. Hence it must thus be the case that its transform causes the V 2 j(p) term 
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to vanish identically too. However, on transforming a(r) = 2d/c 2 r we find that it transforms 
into j(p) = 2d/c 2 p, with insertion of this form for j(p) into V 2 j yielding the non-vanishing 
V 2 (2d/c 2 p) = — (8ird/c 2 )6(p)/p 2 . However, this contradiction is only an apparent one since 
the quantity V 2 (2d/c 2 p) only becomes singular at p = 0, a region where j(p) as given by 
2d/c 2 p becomes too large for a linearized approximation to hold in the first place. The 
apparent non-vanishing of V 2 j is thus a spurious artifact of the linearization in the isotropic 
coordinate system, with j(p) being able to be equal to 2d/c 2 p for all not too small p, and 
with the V 2 j(p) term which appears in the Ricci scalar term in the left-hand side Eq. ( |37|) 
then not coupling to the p(p) source which appears on its right-hand side at all. Thus given 
this pecularity of the isotropic coordinate system, we must solve Eq. (I3"T|) by setting the 
V 2 j term equal to zero in it (while likewise setting the V 4 j term equal to zero in Eq. (13151) ). 
to then reduce Eqs. fl37|) and fl35l) to 



S lv 2 h = P ^ -^h=^, (65) 
6 c 3 c 

to thereby yield attractive gravity after all; with the full solution to the isotropic coordinate 
system Eqs. (1351) - (|37|) in the region exterior to a source of radius p then being given by 

A . 3Mcp 3v , , A . 3Mc ., A . 3Mc 3Mcp 
327To; g 4a 9 p 27r5op 27r5op 32na g 

in complete analog to the standard coordinate system solution given previously. 

Pedagogically, it is instructive to explore this peculiarity of isotropic coordinate system 
linearization a little further. If we return to the full non-linearized theory, we can write 
closed form expressions for the G t t component of the Einstein tensor in both the standard 
(STA) and the isotropic (ISO) coordinate systems, to respectively obtain 

(j* t (STA) = A + J_ _ _L (67) 

and 

V 2 7' 3 V 2 

G'.dSO) = 4 - M. + |L, (68) 

with these two expression being related by the exact coordinate transform given in Eq. 
fl62|) . a transformation which converts the first derivative function G* 4 (STA) into the second 
derivative function G t t (JSO). (In Eq. (16 7p the primes denote derivatives with respect to r 
and in Eq. f[6"81) they denote derivatives with respect to p.) For our purposes here it is more 
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convenient to replace J(p) by the form 



J{p) = [l + K{p)}\ 
as the isotropic coordinate coordinate G* f (ISO) then gets rewritten as 



(69) 



G" t (ISO) 



4 



1 + Kf 



K" + 



2K' 



P 



(70) 



If we now consider an Einstein theory with no source at all, i.e. if we simply set G t t 
equal to zero everywhere, we then find exact solutions to the equation G t t = in the two 



coordinate systems to be of the form 

A{r) 

J(p) = (l 



2C 



1 - 



-) , K(p) = -, 
2pj ' KHJ 2p 



(71) 



(72) 



where C is an unspecified constant and where the two solutions are characterized by the 
same parameter C as they are related by the exact coordinate transform given in Eq. ( 1621) . 
That this form for K(p) really is a solution to G t t (lSO) = is due to the fact that while 
the action of the K" + 2K'/p term on C/2p gives the non-zero —2Cir6(p)/p 2 , the prefactor 
1/(1 + K) 5 vanishes as p 5 near p = to then cancel the delta function. However, suppose 
we now linearize the theory around flat spacetime to obtain 



G^STA; LIN) = — + ^— ^ 



G 4 t (ISO;LIN) = -4 



K" 



2K' 
P 



Setting G* 4 (STA; LIN) equal to zero would still possess a solution for all r, viz. 



A(r) 



2C 

r 



(73) 
(74) 

(75) 



but the function into which it would transform via a coordinate transformation between the 
two coordinates systems, viz. 



9G G 
J(p) = l + T , K(j,) = -, 



(76) 



would not satisfy K" + IK' j p = at p = because of the delta function term that it 
would generate in it. Consequently, even while the solution of Eq. (1751) is a solution to 
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G t t (STA; LIN) = for all r, the function into which it transforms is not a solution to 
G t t (lSO; LIN) = for all p, with the reason for this being that in its linearization G* t (ISO) 
loses its crucial 1/(1 + K) 5 prefactor. Thus, with the coordinate transformation of Eq. 

relating derivatives of different orders, linearization in an isotropic coordinate system 
is not as straightforward as linearization in a standard coordinate system, with it being the 
linearized standard coordinate system which is best suited to accommodate singularities. 

VI. ON THE PHYSICAL SIGNIFICANCE OF THE SCALAR FIELD 

While we have now achieved our primary purpose of showing that the standard 
Schwarzschild phenomenology can in fact be recovered in conformal gravity in the presence 
of a macrosopic scalar field, it is nonetheless instructive to discuss the physical significance 
of the scalar field and its relation first to dynamical mass generation and then to dynamical 
localization of particles in conformally invariant theories. It is also instructive to explore 
the circumstances under which the —S 2 R IX „/V2 term in the action of Eq. (j2J) could in fact 
lead to repulsive gravity. While we have, for simplicity, used a fundamental scalar field for 
the analysis presented in this paper, in the dynamical case the scalar field would be the vac- 
uum expectation value of a fermion composite operator such as the fermion mass generating 
bilinear , ijj(x)i()(x). And whether or not the expectation value of any such composite is to 
be constant or spatially dependent will depend on the state in which the expectation value 
is to be taken. Moreover, with the discussion of dynamical symmetry breaking by fermion 
composites having been formulated primarily in flat spacetime, we will need to adapt the 
formulation to curved spacetime. 

With regard first to the discussion of dynamical symmetry breaking in flat spacetime it- 
self, following the pioneering work of Nambu and Jona-Lasinio [8j based on analogy with the 
BCS theory of superconductivity, we recall that in order to find the ground state of an inter- 
acting fermionic system one looks at states in which all the negative energy fermionic states 
are occupied (for illustrative purposes in the following we shall denote the number of such 
occupied states as N), with all the positive energy fermionic states being unoccupied. Two 
candidate ground states are suggested: the normal one \N) in which the fermion is mass- 
less and the bilinear expectation value (N\ip(x)i(j(x)\N) is zero, and the self-consistently 
determined spontaneously broken superconducting type state \S) in which the fermion is 
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massive and the bilinear expectation value (S\ijj(x)ip(x)\S) is non-zero. Dynamical symme- 
try breaking then occurs in a theory with no fundamental scale whenever the dynamics is 
such that the state \S) has lower energy than the state \N), with the non- vanishing of the 
expectation value (S\ijj(x)ip(x)\S) then introducing a scale dynamically. If such a state \S) 
is to be the vacuum of the theory, then with the ground state needing to be translation in- 
variant, the expectation value (S\~if)(x)if)(x)\S) would then be independent of the coordinate 
x M and thus be constant throughout the spacetime. As such, the lowering of the energy 
from that in the state \N) (viz. (N\H\N) where H is the Hamiltonan) to that in the state 
\S) (viz. (S\H\S)) would yield a cosmological constant term which while unimportant in 
flat spacetime (where only energy differences are observable) becomes central once gravity 
is introduced since gravity couples to energy itself. 

With respect to such a state \S) one can introduce creation and annihilation operators 
for the fermion (generically and a), with the annihilation operators then annihilating the 
state IS") according to a\S) = 0. Given such a Fock space it is very tempting to identify the 
N +1 particle state a^\S) as the lowest lying positive energy fermionic mode, a mode which 
propagates as a translation invariant plane wave above a Dirac sea filled with N negative 
energy plane wave modes. However, such an identification is not necessarily correct since the 
state \S) was determined to be the self-consistent N particle state, with the state (v\S) not 
automatically being a self-consistent N + 1 particle eigenstate of the Hamiltonian. Rather, 
the N + 1 particle state needs to be calculated self-consistently all over again, with it 
beingfound in certain dynamical cases (the double-well potential fundamental scalar field 
case [9, [h]] and a Nambu-Jona-Lasinio type fermion composite model 11|) that the N + 1 
particle state could lower its energy with respect to the state a'\S) by becoming spatially 
dependent, with this spatially dependent state then being the self-consistent N + 1 particle 
eigenstate of the Hamiltonian. Specifically, one constructs a spatially dependent N particle 
coherent state \C) in which the expectation value (C\ip(x)i()(x)\C) is spatially dependent 
(|C) is constructed from \S) via a spatially dependent Bogolubov transform on it), one 
defines creation and annihilation operators (generically and b) with respect to \C) so that 
b\C) = 0, and then constructs the N + 1 particle state w\C). Such an iV+1 particle state is 
not built out of plane waves. Rather, all the N negative energy modes are distorted in the 
vicinity of some center of localization, with the lowest positive energy state then being bound 
to these distorted negative energy states with a wave function which is highly localized and 
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which falls off very fast as we go away from the center of localization. Moreover, the state \C) 
itself need not be an N particle eigenstate at all. Rather, it is the iV+ 1 particle b'\C) which 
is the eigenstate, being stabilized through a cooperative effect between the positive energy 
fermion and the filled negative energy sea in which the localizing of the positive energy 
fermion distorts the wave functions of the fermions in the negative energy sea causing them 
to form a localized potential in which the positive energy fermion is then bound. With 
the spatially independent state \S) being the ground state of the system rather than the 
spatially dependent state \C), there is no violation of translation invariance, with it thus 
being possible to produce localized excited states in a conformal invariant theory without 
needing to break translation invariance. 

While the negative energy modes are distorted in the vicinity of the center of localiza- 
tion in such coherent states, far from the center of localization their wave functions revert 
back to plane waves. Consequently, asymptotically far from the center of localization the 
spatially dependent expectation value (C\ip(x)ip(x)\C) will approach the constant value as- 
sociated with the expectation value (S\ijj(x)ip(x)\S). If we thus represent (C\xjj(x)ip(x)\C) 
by a spatially dependent order parameter S(x) and represent (S\ijj(x)ip(x)\S) by a spatially 
independent order parameter So, we see that S(x) will approach S asymptotically far from 
the center of localization. Moreover, not only will such an S(x) approach So asymptotically, 



in the flat spacetime dynamical symmetry breaking models studied in [11|, [12J the approach 
is exponentially fast (c.f. the dynamical kink type states of 111] where S(x) ~ Sotanh(M;r), 
with S(x) behaving asymptotically as So(l — 0(e~ 2Mx )) where M is the dynamically induced 
fermion mass). Thus while S(x) itself might be large, its derivatives would fall off very fast. 
This then is the localization mechanism for dynamical symmetry breaking in flat spacetime. 

In extending these ideas to curved space, we need to interface them with early universe 
cosmology since mass generating phase transitions occur as the universe cools down, with 
the fermion masses and the fermion condensate S(x) being thought to be associated with 
the electroweak symmetry breaking phase transition which occurs at a temperature Tew 
of the order of 10 15 °K or so. However, even above such temperatures there are already 
dynamical mass scales present such as for instance those associated with the expansion radius 
and scalar curvature of the universe. To generate scales such as these there should thus 
also be a much higher temperature early universe phase transition at some generic grand- 
unified type scale temperature Tqut of the order of 10 28 °K or so. Since this much higher 



20 



temperature phase transition does not generate fermion masses, its associated condensate 
should not Yukawa couple to fermions, with a natural candidate for the needed condensate 
thus being the expectation value of a fermion quadrilinear ip(x)ijj(x)i(j(x)il ! {x) in some 'ur' 
state \U) with expectation value U(x) = (U\'i(j(x)il!(x)'i(j(x)'ijj(x)\U) . The state \U) thus 
has the property that at temperatures above T EW the expectation value {U\?(}(x)il>(x)\U) of 
the fermion bilinear is zero. While these general remarks hold in both standard cosmology 
and conformal gravity cosmology, in the conformal case, just like S(x), U(x) needs to be 
conformally coupled to gravity with the full energy-momentum tensor needing to be traceless. 
In the presence of the two fields U (x) and S(x) we thus replace the conformal matter action 
I M of Eq. © by 

1 1 



dH-g) 112 

1 



2 s^s,, - ^s 2 r^ + as 4 + ijrffrM + rv(x)]v - hSW 

(77) 



2 ,M 12 

with the full T^ u of Eqs. ([7j) and ([8]) being generalized to 

= -\(p + p)U^U v + Pg n + \s^S' u - lg^S' a S. a - IsS'^ + lg^SS' a a 
c 3 6 3 3' 

--S 2 ( - -g^R a Q ) - g' lu \S A + -U'W - -g^U' a U. a - -UU^ V 
6 V 2 / 3 6 y '3 

+ l g ^ UU ;» a - (V - \g^R a ^ - g^rU\ (78) 
and with Eq. (TTOT) being extended to the two scalar field equations 

h 

6 

In the explicit application of conformal gravity to cosmology, we need to take into con- 
sideration both the background cosmology and the fluctuations around it. The background 
cosmology is homogeneous and highly symmetric being either Robert son- Walker or de Sit- 
ter (for simplicity in the following we shall take the background to be de Sitter), while the 
fluctuations around it are inhomogeneous and far less symmetric. (While conformal grav- 
ity has no apparent need for inflation since its Robertson- Walker phase is free of both the 
flatness and the horizon problems [1], a de Sitter cosmology is nonetheless still an allowed 
solution in the conformal case.) The background geometry can thus be modeled by constant 



S 2 R\ = - [p - 3p] + 4AS 4 - (79) 
U 2 R\ = 4tU 4 - UU'^. (80) 
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values for the U(x) and S(x) fields, a background in which the fermions propagate as plane 
waves, while the inhomogeneities can be described by fermions localized in coherent states 
in which the positive energy fermion is localized into some finite region r < R and in which 
the negative energy fermion wave functions only approach plane waves asymptotically far 
from the center of localization. However, since both the fermion bilinear and quadrilinear 
condensates are built out of one and the same set of fermions, the distortion of the fermion 
wave functions in the vicinity of the center of localization will cause both U(x) and S(x) to 
depart locally from their asymptotic background U and So values. Since in principle these 
two departures are not the same, and since the two fields transform independently under 
a conformal transformation, any local conformal transformation which might bring one of 
them to a constant throughout the spacetime would not simultaneously bring the other one 
to a constant form as well. We thus cannot take both U(x) and S(x) to be constant, and so 
in the following we shall in fact allow both of them to acquire a spatial dependence once an 
inhomogeneity is introduced, and shall use these spatial dependences to explicitly monitor 
the modifications to the geometry caused by the introduction of the inhomogeneity. With 
the de Sitter background geometry actually being writable in the static coordinate form 
given in Eq. (1201) . for the case of a single mass source localized in a de Sitter background 
the equations of motion of Eqs. (1221) and (|29|) thus generalize to 
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(81) 

(82) 
(83) 



to thereby define the model. As we see, the great utility of allowing both U(r) and S(r) 
to be spatially dependent is that we can take full advantage of the conformal structure 
of the metric of Eq. (fl9|) to obtain equations which are linear in the metric coefficient 
B(r), regardless in fact of the magnitude of B(r), and thus regardless of whether or not the 
geometry might be close to flat. What prevents finding analytic solutions to these equations 
is only that they are non-linear in the scalar fields. 

With a background de Sitter (or Robertson- Walker) geometry being conformal to flat, in 
such a geometry the Weyl tensor vanishes, and thus also the tensor W^ u introduced in Eq. 
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(J6]). In a de Sitter geometry then, if we set take both U(r) and S(r) to be constant, we find 
that Eqs. flHTT) - fl83|) admit of the exact de Sitter geometry solution 

B(r) = l-kr 2 , Po+Po = 0, -2^ 2 =4Ag 4 + (P °~ 3Po) , - 2kU 2 = 4rE#, (84) 

where the spatially independent po and po denote the background values of p and p. (The 
de Sitter solution is an example of a solution to Eqs. ( 1HTT) - (1831) in which the geometry is 
far from flat even though the equations of motion are linear in B(r).) As such, the solution 
of Eq. fl84|) not only defines the background, in it we see that the explicit removal of the 
conformal factor p 2 (r)/r 2 in going from the metric of Eq. (TT9j) to that of Eq. (|20|) then gives 
us a form of the metric which precisely is de Sitter, viz. 

ds > = _(i _ kr 2 )c 2 dt 2 + " , 2 , + r 2 dn 2 , (85) 

(1 — kr z ) 

in the solution. Use of the metric of Eq. ( 1201) thus precisely dovetails with the use of constant 
background scalar fields. And, moreover, as we see, one should not drop the 4AS'q term as 
had been done in Eq. (TTTT) since this term supports the —(S 2 /6)(R flL/ — \g^ u R a a ) term that 
was retained. Finally, we note that since the Schwarzschild-de Sitter metric can also be 
written in the form of the metric given in Eq. (1201) (with B(r) = 1 — 2/3 /r — kr 2 ), Eqs. (IHTj) 
- ( l83l) are thus ideal for monitoring inhomogeneities in a de Sitter background. 

However, before doing this, we note that since W^ v does vanish in a de Sitter or 
Robertson- Walker background, then, according to Eq. ([5]), in such a background the full 
T^ u of Eq. (1781 must vanish also, with its vanishing then yielding the equation of motion. 



l -[U% + S 2 ] (rt - \g^R a }j =\[{p + v)U^ + pgT\ - g^rUt - g^\S 4 (86) 

when the scalar fields assume their background values. As such, we recognize Eq. flHBI) to 
be in the form of none other than the standard Friedmann cosmological evolution equation, 



save only that the standard attractive Newton constant has been replaced 13j] by a global 
effective repulsive gravitational constant Gglob of the form 

SttCglob 2 

3c 3 " (u$ + s 2 y [ * n 

In conformal gravity then wrong sign terms such as -^R^JYI and -S 2 R^ fl /12 really do 
lead to a gravity which is repulsive, but onglobal rather than local scales. Moreover, this 
was found to be an advantage as it leads [l| to a cosmology with no flatness problem, no 
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horizon problem, no universe age problem, to a cosmology which is naturally accelerating 
without fine-tuning, and to a solution to the cosmological constant problem; with it not 
being the cosmological constant which is quenched down from the large value that would 
be associated with a large Uq and So, but rather it is the amount (viz. Gglob) by which 
it gravitates which is quenched instead, as the very same mechanism which makes Uq and 
Sq large simultaneously makes Gglob small, with the conformal theory then being able to 
naturally fit [1] the accelerating universe supernovae data without any fine-tuning at all 
despite using a cosmological constant which is as large as particle physics suggests. 

With mass sources only being able to localize in the presence of inhomogeneities, the local 
gravity produced by their presence will be associated with departures of the scalar fields 
from their constant background values and with departures of the Weyl tensor from its zero 
background value. In such a situation local gravitational departures from the background 
geometry will be controlled by the sign of a g , with there then being an induced local Gloc 
whose sign is completely decoupled from that of the global Gglob, with the — £7 2 -R M M /12 and 
—S 2 R fJ, /1 /12 terms thus not leading to repulsive local gravity. To see exactly how it works 
in the conformal theory we thus need to look at the non- leading terms in U(r) and S(r), 
terms which will explicitly be spatially dependent. 

While we had noted above that in the flat space case departures from constant scalar field 
are exponential, an explicit study of Einstein gravity coupled to a scalar field has found [14J] 
that in the gravitational case such departures are only power behaved. Specifically it was 
found that an Einstein gravity theory coupled to the kinetic energy of a minimally coupled 
massless scalar field, viz. one based on the action 



\ n + \d"Sd,S 



(88) 



and equations of motion 
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or 



admitted the exact exterior isotropic coordinate system solution 
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where the parameter d is given by 

d = |l6p 2 - ^— ) (91) 
24 



and where C, p and So are appropriate integration constants. As such this solution itself 
generalizes the point source solution to the same theory found even earlier by Yilmaz [3] in 
which the integration constant p was zero, with the non-vanishing of p in Eq. (M?j) enabling 
the solution to describe an extended source with a radius of order p rather than a point 
source. On setting the parameter d equal to 2MG/c 2 , and converting the solution to the 
standard coordinate system of Eq. (|39|) . at large values of the radial coordinate r the metric 
coefficients are found to behave as 

2MG AMGp 2 M 3 G 3 

B[r) —> 1 1 , 

c 2 r 3c 2 r 3 3c 6 r 3 ' 

1 1 2MG , 4p 2 M 2 G 2 AMGp 2 M 3 G 3 



to thus yield power law corrections to the Schwarzschild metric. Moreover, it was further 

n 

shown in |14j that if the action of Eq. (188j) were to be augmented by a double-well Higgs 
potential of the form V(S) = XS 4 — p 2 S 2 /2 with its familiar minimum at S = /i/(4A) 1 / 2 , the 
corrections to the now Schwarzschild-de Sitter metric would continue to be power behaved. 
With such corrections being only power behaved rather than exponential, it was further 
noted in ljj] that if a source such as the sun actually did get its mass dynamically via a 
macroscopic Higgs scalar field, the Higgs mechanism could actually be tested gravitationally 
by monitoring these power law corrections to the Schwarzschild geometry, as the corrections 
associated with the expansion of Eq. fl92|) could potentially be quite significant if the pa- 
rameter p is substantially larger than the Schwarzschild radius of the sun, something that 
would indeed be the case if p is of order the radius of the sun. 



VII. POWER-BEHAVED SOLUTION TO THE FULL CONFORMAL THEORY 

Given the above remarks, we shall thus seek power-behaved corrections to the conformal 
gravity background solution of Eq. (1541) . and since we are in the fourth order theory we 
shall need to incorporate the linear potential term of the solution of Eq. (1231) as well. We 
thus look to see whether Eqs. (15TT) - (183"]) can support a power-behaved exterior solution of 
the form 

B(r > R) = -kr 2 + <yr + w- ^ + ^ + ..., 

Q O O O O 

S(r> R) = S + — + — + — + — + — + ..., 

fy fy>£ iyO ny L ± <y<j 
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U(r > R) = U 



r -. 1 2 u 3 r, u ; 



ry ,y2 ry^ i^>4 



+ . . . 



(93) 



in the event that there is some localized source within the r < R region, a source whose 
very presence will also lead to exterior region corrections to the background global p(r) and 
p(r) (by distorting the plane wave modes of the negative energy fermion sea). As we shall 
see, these corrections to the perfect fluid energy density and pressure will be power behaved 
too, with the full energy density and pressure behaving as 



p(r) = 9(R - r)p LOC (r) + 9(r - R)p GLOB (r) 

ry ry 2 ry 3 ry^. i^>5 

p(r) = 9{R - r)p hOC (r) + 9(r - -R)pcLOB(r) 



6(R-r)p LOC (r) + 9(r-R) 



9{R - r)p hOC {r) + 9{r - R) 



. Pi , P2 . P3 . Pi . P5 . 

Po + — + — + ^ + — + ^ + 

ry ry _■ ryO ty~) 



(94) 



where pLoc( r ) and Ploc(^) are the localized r < R contributions to the energy density and 
pressure due to the presence of the localized positive energy fermion. 

To see how the theory is able to support power-behaved corrections we need to insert the 
expansions of Eqs. (l9~3l) and ( 1941) into the relevant quantities which appear in Eqs. (1511) - 
to obtain first for Eq. (1521 (and analogously for Eq. 
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4:kS S 3 + 2 1 S S 2 + 2(W ^ SoSi 
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8kS S A - 2kS 1 S 3 



+ 5 7 S S 3 + 7S1 S 2 + 2{W n ^ S S 2 + 



-S? + 2wS S 2 - 2/3S Si 
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(w - 1) oL > 

-|- 4WOQ02 — ^A- 

UkS S 5 - QkS^i - 2kS 2 S 3 + W-fS S 4 + 4rfSiS 3 + 7 S 2 2 + 



2(w 



-SnS; 



CM3 



2(w 



4AS 4 = 4AS# + 



16A5 3 5i 



-StS 2 + 6wS S 3 + 2wS 1 S 2 - 8f3S S 2 - 2(3Sl + 25S , S'i 



(95) 



1 

+ ^7 



1 



+- 



+- 



16\S$S 2 + 24\S$Sl 
16\S*S 3 + 48AS 2 S 1 S 2 + 16XS Sf 
16AS 3 S 4 + 24AS 2 S 2 + 48AS 2 S 1 S 3 + A8XS S 2 1 S 2 + 4\S* 
16\S*S 5 + 48 AS 2 ^4 + 48\S*S 2 S 3 + 48AS SiS 2 2 + 48A5 S 2 S 3 + 16XSfS 2 



(96) 
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o o i (Pi - 3 Pl) i (P2 - 3p 2 ) , (p 3 ~ 3p 3 ) , (P4 - 3p 4 ) , (P5 ~ 3p 5 ) 

p - 3p = p — 3p H 1 o 1 s 1 i 1 is ' 

(97) 

as evaluated in the r > R region. In the above we have carried the expansions as far as 
needed in order to exhibit at which particular level each of the explicitly displayed terms in 
Eqs. ( 1931) and ( 1941) first appears. With Eqs. ( 1951) - ( 1971) all displaying the same power series 
structure, we see that Eqs. ( 182]) and ( 1831) will indeed support a power series expansion order 
by order. The zeroth order term is already satisfied by the background solution given in Eq. 
( 184|) . with the next few orders being maintained by 

-2kS S 1+1 S 2 = 16Aff 3 S 1+ (pi ~ 3pi) , 
-2fcC/ £/i + 7^ Q 2 = lerf/g 3 ^, 
-2A;5 5 2 + 7^1 + ^ W ~^ S 2 = 16XS 3 S 2 + 2AXS 2 S 2 1 + (P2 ~ 3P2) , 

o C 

— 2fcf7 ^2 + 7 Wi + (W 7 ^ U 2 = 16tU 3 U 2 + 2AtU 2 U 2 , 

o 

-4kS S 3 + 2yS S 2 + -^-p—^-SoSt = 16XS 3 S 3 + 48XS 2 S 1 S 2 + 16XS S 3 l + ( ^^, 

-4kU U 3 + 2 7 £W 2 + 2(W ~ ^ CWi = 16rf/ 3 t/ 3 + 48rf/ 2 [/ 1 ^2 + 16^0^. (98) 

The pure C/(r) sector solution can be simplified to 

773 

k = -2rU 2 , 7 = UtUqUu w-1 = 36tU U 2 + 36rU 2 , U 3 = -|, 

/3 = 2C/ C/ 2 + 72rU U 2 U 2 + 72rU 2 U 2 + 24rU?, (99) 

with the S(r) sector solution then being satisfied by adjusting the p — 3p power series terms 
accordingly. 

For the treatment of Eq. (1HT1) . we note that in it there is a factor of B(r) in the de- 
nominator. However, for a de Sitter geometry the coefficient B(r) can actually vanish, 
with the radial distance r then being bounded by some maximum cut-off value A. (With 
the large r limit of the candidate B(r > R) metric given in Eq. (1931) being of the form 
B(r > r) — > w + 7r — kr 2 , positivity of the parameters w, 7 and k for instance would entail 
a cut-off at A w (y/2k)[l + (1 + to/7 2 ) 1 / 2 ].) For the function f(r) given in Eq. ( 18 ID to be 
well-defined, the quantity p + p would need to vanish at r = A. Now this quantity actually 
does vanish at r = A in the de Sitter background, since as noted in Eq. ([HID , Po +Po actually 
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vanishes everywhere in the de Sitter background. We thus now need the full p + p to vanish 
at r = A also. Since the p + p and p — 3p combinations are linearly independent, we are thus 
free to adjust p + p without needing to affect the treatment of Eq. (1521) which we have just 
given. To achieve the needed vanishing of p + p at r = A we must set it equal to (A-r)D(r) 
where D(r) is a function of r which regular at r = A. (In terms of a p + p expansion of 
the form J2(Pn +Pn)/r n , this is equivalent to setting J2(Pn +Pn)/A n = 0.) With B(r) also 
possessing a simple zero at r = A, we can expand the {p + p)/ B ratio as a power series in r 
in the r > R region, to obtain 

3[p(r >R)+p(r > R)] = ^ + R, + B* + R* + R, + R, 
cB{r > R) r r 2 r 3 r 4 r 5 

With the existence of the cut-off, the solution to Eq. ( 1511) now takes the form 

S(r) = -- f dr'r' 2 f(r') - — f dr'r ,A f(r') - - dr'r' 3 f(r') - — dr'r' J Or'), (101) 
2 Jo 6r Jo 2 J r 6 JV 

as now integrated up to A. 

For the candidate expansion of Eq. (1931) the scalar field term which appears in the source 
function f{r) in Eq. (15T|) is given by 

^SqSi QSqS2 125'o5'3 (205 , q5 , 4 -I- 25'x>S'3 — 25*2^ 



SS" - 2S' Z = -^p + + — + 



y*3 ^4 jyi5 



| (SOgogs + 63^4 - 6 5*25*3) 

together with an analogous expression for the L r (r)-dependent term. With the expansion of 
Eq. (I102p only beginning in order 1/r 3 , we see that there is nothing available to cancel the 
Ro, Ri/r and R2/r 2 terms which appear in Eq. ( 11001) . With the insertion of these three 
particular terms leading to contributions to 5(r > R) which respectively behave as r 4 , r 3 
and r 2 logr, to recover the candidate expansion of Eq. ( 1931) . we must thus set Ro = 0, R\ = 0, 
and i?2 = 0. Additionally, we note that each one of the 1/r 3 , 1/r 4 and 1/r 5 terms which 
appear in Eqs. (11001) and H102[) would also generate a term in B(r) which involves logr. To 
recover our candidate form for B(r > R) we must thus cancel all such terms as well, and 
must thus set 

- J R 3 + 25o5 1 + 2f/ f/i = 0, - J R 4 + 65'o5 2 + 6f/o^2 = 0, -i?5 + 125 5 3 + 12f/o£/3 = 0. (103) 

While we need to cancel the contributions of the 1/r 3 , 1/r 4 and 1/r 5 terms, for the higher 
powers everything is well behaved. And not only that, the insertion of such higher powers 
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into Eq. (11 Oil) will generate back for us the candidate power series for B(r) precisely as 
given in Eq. fl93l) . with the insertion into Eq. (110 II) of the source term in the generic form 
f(r >R) = E 6 °° f n /r n yielding 



B(r > R) = J2 fn 



n=6 
+ 



1 1 r 2 

+ 777 7T7 7 + 



2(n-3)R n - 3 6r(n-5)R n - 5 2(n - A)A n ~ i 6(n - 2)A™~ 2 
1 



r«- 4 (ra - 2)(n - 3)(n - 4)(n - 5). 
just as desired. Finally to get the full solution, on setting 



(104) 



3 

hoc = ~ 4a c B <KphOC + PhOC ^ ( 105 ) 
we need to augment the globally generated Eq. (11041) with the additional locally generated 

B(r >R) = ~ [ R dr>r> 2 f LOC (r') - 1 [* dr'r' A f LOC (r'), (106) 
2 jo or jo 

to thus yield a net 1/r term in B(r > i?) of the form 

With regard to Eq. (I107p . we note that even though it is difficult to assess the relative 
importance of the various terms which appear in it, with the entire source term of Eq. (1HT1) 
being multiplied by an overall l/a g factor, we can nonetheless conclude that there will be 
some choice for the sign of a g for which the net effect in Eq. (11071) will be attractive just as 
desired, with the signs of the various terms which appear in Eqs. (|98|) and (|99|) adjusting 
accordingly. As we thus see, an attractive 1/r potential can naturally arise in conformal 
gravity even when long range macrosopic scalar fields are present. 



VIII. MODIFICATION OF TRAJECTORIES DUE TO THE SCALAR FIELD 

With the scalar field modifications to the exterior metric associated with Eqs. (11041) and 
(11061) beginning in order 1/r 2 , a precise enough monitoring of the geometry in the vicinity 
of the sun could thus reveal the presence of any macroscopic scalar field that the sun might 
possess. For any monitoring which involves photons, we note that since massless particles do 
not couple to the mass generating S(r), for light rays one can continue to use the standard 
massless test particle geodesies, as evaluated now in the metric associated with Eqs. ( 11041) 
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and (I106p . For massive particles however, there is an additional effect, since their very 
coupling to the spatially dependent S(r) itself also modifies their motion. To study this 
effect in detail we would need to eikonalize the Dirac wave equation given in Eq. ([3]) to find 
the trajectories of its short wavelength ray solutions. To get a sense as to what such eikonal 
trajectories might look like, we can instead vary the test particle action 



-h J drS(x), (108) 



as this action not only reduces to the conventionally used massive test particle action when 
S(x) is constant, but for varying S(x) it is actually fully conformal invariant, since the e~ a ^ 
change in the scalar field is compensated by an accompanying e a<yX ^ change in the proper 
time. Variation of the action of Eq. (I108P with respect to the coordinates of the test particle 
thus yields conformal invariant trajectories, with their specific form being given as 17] 



, „ ( d 2 x x , dx^ dx v \ ( , n dx x dx f3 \ 

ks (TP" + r -lFlfr ) = [*" + H7H7)- < 109 > 

With the Yukawa coupling constant h canceling in Eq. (11091) . and with all fundamental 
fermions getting their masses from one and the same SU(2) x U(l) Higgs doublet (with 
masses that only differ because each fermion has its own associated h), we see that even 
with a spatially varying S(r), all particles will nonetheless still fall in a gravitational field 
with trajectories which are independent of their masses. 

To determine the specific trajectories which are involved in the static, spherically sym- 
metric source case of interest, we note that for the general standard coordinate metric given 
in Eq. (1391) . the four equations of motion contained in Eq. (I109j) take the form 

dH B' dt dr S' dt dr 
dr 2 B dr dr S dr dr ' 

d 2 r A' fdrV r ( d6\ 2 rsin 2 fl ( dchV B' ( dtV S' S' ( drV 



dr 2 2A\dr) A \dr J A \dr J 2A \dr J AS S \dr J 

d 2 6 2 d6dr _ sinQcose fd£\ 2 _ S'dOdr 

dr 2 r dr dr \ dr J S dr dr' 

d?4> 2 d(f> dr ^ ^ ^dcj) d6 S' d<p dr (110) 

dr 2 r dr dr dr dr S dr dr 

As we thus see, even in the presence of a spatially dependent scalar field, the equations of 
motion still admit of trajectories with fixed 9 = tt/2. In such trajectories we find that the 
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three other equations of motion admit of exact first integrals 

BS^- = C, 

(IT 

, ,(dr\ 2 K 2 C 2 
S+AS {Tr) + ^~T = D ' 

r 2 sf T = K, (111) 

where C, D and K are integration constants. Finally, on eliminating the dependence on dr, 
we find that the trajectories can be written as 

c2 AC 2 fdrV K C 2 

s2 + AK^(dr_Y + ^_C^ = 

r 4 \d*f>) r 2 B 2 ' 

r 2 dd K . . 

51 = o- (112) 

When S is constant these relations all reduce to the standard massive particle trajectories, 
with the modifications when S is spatially varying then providing a possible window on 
symmetry breaking physics as well as constraints on the magnitudes of such modifications. 
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